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AN APPROACH TO THE DETERMINATION 
OF AIRCRAFT HANDLING QUALITIES BY USING 
PILOT TRANSFER FUNCTIONS 

By James J. Adams and Howard G. Hatch, Jr. 

Langley Research Center 

SUMMARY 

It is shown that a correlation exists between longitudinal closed-loop characteristics 
of a pilot-aircraft system (determined by using analytical expressions for pilot response 
along with the analytical expression for the aircraft response) and pilot ratings obtained 
in many previous flight and simulation studies. Two levels of preferred pilot response 
used in the computations are (1) a static gain and a second-order lag function with a lag 
time constant of 0.2 second; and (2) a static gain, a lead time constant of 1 second, and a 
lag time constant of 0.2 second. 

If a system response with a pitch-angle time constant of no more than 2.6 seconds 
and with a stable oscillatory angle -of -attack mode of motion having a period of no more 
than 2.5 seconds can be achieved with the first-level pilot model, the pilot rating for the 
vehicle is satisfactory. Further, if an altitude response with a stable oscillatory mode 
of motion with a period of no more than 5 seconds can be achieved, the vehicle is rated 
satisfactory. If the second-level pilot model is required to achieve these system response 
characteristics, the vehicle is rated acceptable-unsatisfactory. 

INTRODUCTION 

It has been shown in references 1 and 2 that models of human response together 
with models or equations of motion of the vehicle can be used to predict system perfor- 
mance of manually controlled vehicles. It would be beneficial if these pilot models could 
also be used to predict pilot acceptance, or rating, of future vehicles. This thought has 
been discussed previously in reference 3. Since a large amount of previously obtained 
data is available on pilot ratings of aircraft as functions of aircraft characteristics, 
criteria for predicting pilot ratings are established in this paper by showing the correla- 
tion between calculated characteristics of the pilot-vehicle system (determined by using 
pilot models) and the previously obtained pilot ratings of aircraft. The criteria that are 
developed involve both specified levels of preferred pilot response characteristics and 



specified closed-loop characteristics of the pilot -vehicle system. The emphasis in this 
paper is on pilot ratings for aircraft longitudinal control, but the procedure should be 
applicable to other control situations. 


SYMBOLS 


Values are given in both SI and U.S. Customary Units. The measurements and 
calculations were made in U.S. Customary Units. 
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c mean aerodynamic chord, m (ft) 

e displayed error, rad 

h altitude, m (ft) 

h c ,h g altitude command and error, respectively, m (ft) 

Iy moment of inertia, kg-m^ (slug-ft^) 

K general gain 
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pilot-model static gain 

outer-loop pilot-model static gain 

inner-loop pilot-model static gain 

normalized lift-force derivative, C T , per sec 

2m 


normalized control-lift derivative, C T , per sec 

’ 2m L 6 e ’ * 


normalized damping moment in pitch due to angular pitch rate, 
per sec 


normalized pitching-moment derivative with respect to angle of 
pV S( ? C m , per sec 2 


2I y ' m «’ 


normalized pitching-moment derivative with respect to elevator 


pV Sl - C m , per sec 2 
2 I y 6 e 

mass, kg (slugs) 

closed-loop period of altitude mode of motion, sec 
closed-loop period of angle -of -attack mode of motion, sec 
angular pitch rate, rad/ sec 
surface area of wing, m 2 (ft 2 ) 

Laplace operator, per sec 
closed-loop time constant, sec 
pilot-model lag time constant, sec 
pilot-model lead time constant, sec 


pVSc 2 
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attack, 


deflection, 
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open-loop aircraft lead time constant, sec 


time, sec 

velocity, m/sec (ft/sec) 
displacement error, m (ft) 

angle of attack, rad (except where otherwise noted) 

flight-path angle, rad 

control deflection, rad 

elevator control deflection, rad 

damping ratio 

damping ratio of altitude mode of motion 
damping ratio of angle -of -attack mode of motion 
pitch angle, rad (except where otherwise noted) 
pitch-angle command and error, respectively, rad 
air density, kg/m** (slugs/ft 3 ) 
closed- loop control frequency, rad/sec 

control-actuator undamped natural frequency, rad/sec (or cps) 
closed-loop altitude frequency, rad/sec 
open-loop aircraft undamped natural frequency, rad/ sec 
closed-loop angle-of-attack frequency, rad/sec 



Abbreviation: 

PR pilot rating 

Single and double dots over symbols denote first and second derivatives, respec- 
tively, with respect to time. 


PILOT MODELS 
Single -Loop Control Tasks 

For single-loop control tasks, the pilot model used herein consists of the following 
elements: a static gain and a lead time constant, which constitute a mentally commanded 
control deflection; and a second-order, critically damped lag function, which represents 
the dynamic response of the arm in executing the command for the control deflection. 

In transfer-function form the model is 

6 K l(l + T 2 S ) 

e (1+TjS) 2 

Typical measured pilot transfer functions showing the response of the pilots when con- 
trolling three different plants are shown in table I. Two examples are given for each of 
the plant dynamics to illustrate the range of coefficients that occur. 

TABLE I 


MEASURED PILOT RESPONSE 


Dynamics 

Pilot transfer function, 5/e 

2 

3.26(1 + Os) . 0.7(1 + 0.67s) 

s 

(1 + 0.25s) 2 ’ (1 + 0.33s) 2 

10 

0.5(1 + 0.3s). 0.36(1 + 0.91s) 

s(s + 1) 

(1 + 0.1s) 2 ’ (1 + 0.18s) 2 

10 

0.3(1 + 0.47s). 0.29(1 + Is) 

s 2 (s + 1) 

(1 + 0.033s) 2 ’ (1+ 0.05s) 2 


Average values derived from all the appropriate tests reported in references 4, 
5, and 6 for the lead and lag time constants are presented in table n. 
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TABLE H 


AVERAGE VALUES FOR LEAD AND LAG TIME CONSTANTS 
OF PILOT RESPONSE 


Dynamics 

Lead time constant, 
T 2 , sec 

Lag time constant, 
T 1 , sec 

2 

0.26 

0.22 

s 



10 

.92 

.17 

s(s + 1) 


10 

.97 

.044 

s 2 (s + 1) 


It can be seen from these tables that when a pilot is controlling a very easy-to-handle 
plant, such as a vehicle with pure-rate dynamics 2/s, the lag time constant of the pilot’s 
response Tj is as large as 0.33 second, with 0.22 second being the average value; the 
pilot-model lead time constant T 2 is often zero, with the average value being 0.26 sec- 
ond. However, it seems that a pilot will readily add more lead because sometimes a 
substantial amount is measured for dynamics of 2/s. If the lag of the plant is increased, 


as for the plant with dynamics of 


10 


s(s + 1) 


, lead is always present, and the average value 


of Tp increases to 0.92 second. When more compensation is required on the part of 
the pilot to control the system because the plant lag is increased further, the pilot-model 
lag time constant is reduced. For the type of plant which a human operator is just barely 

/ s. 

10 


able to control in a stable manner I that is, a plant with dynamics of 


, the pilot - 


\ s 2 (s + 1)/ 

model lag time constant is reduced to 0.044 second. It appears from these data that 
appropriate steps in pilot adaptation to plant difficulty are first to add lead and, second, 
to reduce lag. 


In one experiment reported in reference 4, pilot response was measured while con- 
trolling an aircraft-pitch-response type of plant. The plant dynamics were 


0 = 10 
6 s(s 2 + 3s + 10) 

and the measured pilot transfer function was 


(2) 


6 




( 3 ) 


6 = 0.86(1 + 0.71s) 
e (1 + 0.14s) 2 

It can be seen from this pilot transfer function that the lag and lead time constants are 
within the ranges noted in the previous paragraph. It is assumed, therefore, that the 
maximum and minimum average values for the time constants given in table II also apply 
to aircraft plant dynamics. 

The human pilot provides a static gain Kj, which is always adjusted to provide 
desired closed-loop system response characteristics. That these system characteristics 
can be identified and defined is shown by the results obtained in reference 4. In refer- 
ence 4 a wide variety of plants, control gains, display gains, and subjects were tested. 

In spite of this great variety of system elements, the dominant closed-loop system 
characteristic (the mode of motion which provides the largest contribution to the time 
response of the system) always fell within a very restricted range of period or time con- 
stant. For the type of plants considered in reference 4, that is, plants of the form K/s 
and K/s^, which have no static stability of their own, the dominant system character- 
istic was an oscillatory mode of motion. The frequency of this oscillation always fell 
within the range of from 2 to 4 rad/sec (that is, the period was around 2 seconds), and 
the oscillation was always stable. From these results, it is concluded that the pilot uses 
the dominant system characteristics as the criteria to set his static gain. 

The dominant system characteristic is not always an oscillatory mode of motion. 
With the aircraft-pitch-response type of plant, the dominant closed-loop system charac- 
teristic was a first-order mode of motion. (In the case represented by eqs. (2) and (3), 
the time constant of this first-order response was 1.7 seconds.) There was also an 
oscillatory mode of motion included in the system response in this case, and this oscilla- 
tory mode was stable. 

From these examples, it is concluded that the system response which the pilot con- 
cerns himself with can be of two different types. If the dominant system characteristic 
is an oscillatory mode of motion, the pilot is concerned with the frequency and damping of 
this mode of motion. If the dominant system characteristic is a first-order mode of 
motion, the pilot is concerned with the time constant of this mode and the stability of the 
other oscillatory modes of the response. In either case the system time characteristic 
(period or time constant) that the pilot strives to achieve is around 2 seconds. 

Multiloop Control Tasks 

If the control task is a multiloop problem, where there is a second variable (such 
as altitude or horizontal displacement) which is controlled by the manipulation of an inner- 
loop variable (such as pitch angle), it is necessary to specify an outer-loop pilot model in 
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series with the inner-loop pilot model. A block diagram of such a system is shown in 
figure 1. In this type of system the single-loop results discussed in the previous section 
are applied to the inner loop. 



Figure 1.- Block diagram of a typical multiloop system. 


Reference 6 shows that in the multiloop control task of lunar module hover, where 
the outer-loop variable is horizontal translation and the inner-loop variable is pitch 
angle, the series arrangement of the pilot model does give a true representation of human 
pilot response. In this task the outer-loop plant dynamics were 


fe _ 5J36 

e s 2 


and the outer-loop pilot model was 

9 c _ 0.009(1 + 9.2s) 
x e (1 + 0.1s) 2 

It was also shown in reference 6 that the lag function included in the outer-loop pilot 
model is not really required to properly represent the pilot. The outer-loop variable 
(horizontal translation) is controlled at such a low frequency that a lag of the magnitude 
of 0.1 second does not influence the closed-loop response. 

Reference 7 has shown similar results for helicopter hover. Again, a series 
arrangement was used to represent the pilot; the required outer-loop pilot model was 


Q 

~^= 0.027(1 + 0.4s)e-°- 08s 
x e 


( 6 ) 


where a small time delay is used instead of a lag function. 

Reference 8 used the series arrangement to represent the pilot in the case of alti- 
tude control for an airplane; the required outer-loop pilot model was 

e c _ 0.4 

h e 1 + 0.1s 1 

In this case, no lead was required and again a very short lag time constant was included. 
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In all three of these references (refs. 6, 7, and 8), the outer-loop variable is con- 
trolled at a lower frequency than is the inner-loop variable. Also, when the outer loop 
is closed around the inner loop, the resulting modes of motion couple with each other, 
and restrictions on the permissible outer- and inner-loop static gains can result. 

Pilot Response Levels 

In the single-loop and multiloop experiments just described, different plant dynamics 
were used to provide different levels of control difficulty, and the pilot response was 
determined. On the basis of these results, different levels of preferred pilot response 
are formulated, and it is hypothesized that these levels will correspond to different levels 
of pilot ratings. The thought that a pilot would prefer to act as a simple amplifier was 
first discussed in reference 9. In the present paper this idea is expanded to cover the 
whole range of pilot opinion. For single-loop, or inner-loop, control, the levels of pre- 
ferred pilot response are as follows: 

(1) The pilot would prefer to operate as a simple amplifier with a lag time constant 
of 0.2 second. This is a very simple response, achieved by the pilot in a relatively 
leisurely manner. It is hypothesized that if suitable system response characteristics 
can be achieved with this pilot response, the system will be judged satisfactory. 

(2) If moderate compensation is required on the part of the pilot in order to stabilize 
the system the compensation takes the form of a lead with a time constant of 1 second, 
with the lag time constant maintained at 0.2 second. In order to supply this lead the 
pilot must perform the additional mental task of differentiating the displayed information. 
It is hypothesized that if suitable system characteristics can be achieved with this pilot 
response, the system will be judged acceptable-unsatisfactory because of the added 
required task of differentiating. 

(3) If further compensation is required, it is supplied by reducing the lag time con- 
stant, with values as low as 0.05 second being obtainable. To achieve this response the 
pilot must not only perform the additional mental task of differentiating the displayed 
signal but must also increase his physical activity to obtain the faster control deflection. 

It is hypothesized that if this compensation is required to stabilize the system, the system 
will be judged unacceptable. 

For outer-loop control, two levels of preferred response are proposed. They are 
as follows: 

(1) The pilot would prefer to operate as a simple amplifier. It does not appear to 
be necessary to specify any particular value for lag time constant in this case, because 
the expected values of this constant will not influence system response; therefore, a value 
of zero is used. If this type of response is all that is required for the outer loop, no 
degradation of pilot rating will result. 
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(2) If compensation is required, it takes the form of lead. Apparently, outer-loop 
lead time constants can be quite large. If lead is required in the outer-loop pilot model, 
some, as yet undefined, degradation in rating will occur. 

For all levels of pilot response, the model static gain is adjusted to provide the 
laigest closed-loop system frequency (the shortest time characteristic) possible while 
maintaining stability. It is expected that for single-loop or inner-loop variables, a time 
characteristic of around 2 seconds will constitute a desirable system response and that 
for outer-loop variables, the time characteristic will be longer than 2 seconds. It will 
be necessary to determine a suitable value for the time characteristic of the outer-loop 
variable. 

Both single-loop and multiloop systems are examined in this study because both 
types of system are represented in typical flight tasks. The single-loop, attitude control 
task corresponds to long-range tracking tasks. The multiloop, altitude control task cor- 
responds to landing -approach tasks and also terrain-following and formation-flying tasks. 

PILOT -RATING DATA 

A large amount of data on pilot rating as a function of various aircraft parameters 
has been accumulated in the past. Perhaps the most fundamental relation is that estab- 
lished between pilot ratings and aircraft longitudinal short-period characteristics 
and 2£w n . The boundary dividing pilot ratings of satisfactory (PR = 1 to 3.5) from 
acceptable-unsatisfactory (PR = 3.5 to 6.5) and the boundary dividing ratings of acceptable- 
unsatisfactory from unacceptable (PR = 6.5 to 9) are shown in figure 2, where i S 

plotted as a function of 2£w n . (Note that the term acceptable -unsatisfactory is shortened 



Figure 2.- Pilot ratings as functions of aircraft longitudinal short-period 

O 

characteristics tc^ and 2 ^ 0 ^, of 1^, and of cu^. 
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and just the word acceptable is used to denote this rating.) These boundaries were taken 
from reference 10 for a landing condition with good (front-side) drag characteristics - 
that is, with drag that increases with increases in speed. 

Pilot rating has also been shown to be a function of 


M, L - L,. ]V1 
a a 

e e 


6,2 


M 


6 , 


( 8 ) 


Since Lg is usually small, 
e 


L 6 , 2 


a 


The relationship of pilot rating to Tg g has been defined in several ways. In refer- 
ence 11, where the flight tasks include both altitude and attitude control, it is defined 
with the use of the parameter Tg 2 a) n’ use the parameter T g causes the data 

to coalesce to the curve shown in the second part of figure 2. The data of reference 11 
indicate that the best pilot ratings are obtained when Tg 2 w n = ^ and that the ratings 
become less favorable for higher or lower values of T 0 2 a> n . However, a question as 
to whether the ratings do really become less favorable for the high values of Tg 2 u> n 
has been raised in other investigations (e.g., ref. 3). 

Pilot rating has been shown to be a function of control-actuator dynamics also. 
Reference 12, which is a simulator study, shows that when second-order elevator actuators 
are inserted into an aircraft longitudinal control system, pilot ratings can be affected. 

For example, as is shown in the third part of figure 2 (taken from ref. 12), if the natural 
frequency of the actuator u>f is less than 3 cps (20 rad/ sec) the rating of the otherwise 
good aircraft is seriously degraded. 

The purpose of the present study is to show that characteristics of the pilot-aircraft 
system, obtained by using the various levels of preferred pilot response defined in the 
preceding section, correlate with the pilot ratings as functions of c and 2£o> n , of 
Tg 2 > and of actuator dynamics. 

AIRCRAFT AND SYSTEMS REPRESENTATION 


The aircraft is represented by the two-degree-of -freedom longitudinal equations 
of motion 


a - 6 = -L a 
a 


( 9 ) 
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and 


9 = M q 9 + M a a + M g fi e (10) 

and by the relationship for altitude 

h = V(0 - a) (11) 

Two different values of L a were used in the investigation - a typical low-speed 
value L^, = 0.585 per sec and a high-speed value L a = 1.3 per sec. The value of 
Mg was arbitrarily assumed to be 1.0 per sec^. The values of Mq and M a were 

e o 

adjusted to provide a wide range of aircraft short-period characteristics to* and 2£u> n . 
An example of the performance of the pilot-aircraft system is given in appendix A. 

ANALYSIS 

Method 

Under the hypothesis of this paper, a given aircraft in a given flight condition is 
rated satisfactory if suitable closed-loop system response characteristics can be 
achieved by using a pilot model with zero lead and an inner-loop (or single-loop) lag 
time constant of 0.2 second. No restrictions are placed on pilot-model static gain, but 
the limitations on suitable system response characteristics must be defined and met. 

One of the limitations used is that all system oscillatory modes must be at least neutrally 
stable. This minimum damping requirement is used because of the indication in previous 
investigations on pilot response that a human pilot will strive to achieve certain closed- 
loop system frequencies even if system damping must be sacrificed. Determination of 
these desired frequencies (or system time characteristics) is needed to complete the 
definition of suitable system response characteristics. For this purpose an exploratory 
analog computer study was made to determine the response of the pilot-aircraft system 

O 

for a wide variety of combinations of to* and 2£cu n , with L a set equal to 0.585 per 

sec. In this process a trial-and-error adjustment was made of the pilot-model static 

gains to minimize the dominant-mode time characteristics while maintaining at least 

neutrally stable system response. The time characteristics thus determined for values 
o 

of w n and 2£w n in the vicinity of the experimental boundaries between pilot ratings 
of satisfactory and acceptable (that is, at PR = 3.5) should then define the minimum 
suitable system response. 

First, the altitude response obtained by using the multiloop system shown in fig- 
ure 3 was determined. With this system there were two dominant oscillatory modes 
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Pilot model 



Figure 3. - System tlock diagram. 


of motion — an altitude mode of motion which had a relatively long period and an angle - 
of-attack mode of motion. There was also a high-frequency, well damped control mode 
of motion which did not contribute in any significant way to the closed-loop time response 
of the system. 

O 

When aircraft with high values of cu* were considered, the determination of the 
best response was very simple and straightforward because the outer loop did not couple 
with the inner loop in any degrading manner and because the frequency of the inner-loop 
angle -of -attack mode of motion was always higher than was believed necessary to satisfy 
a pilot’s desire for adequate performance. That is, the period, or time characteristic, 
of the angle -of -attack mode of motion was less than 2 seconds. Therefore, the inner- 
loop static gain was adjusted to be as high as it could be without making the angle-of- 
attack mode of motion unstable. In this way the maximum amount of compensation was 
provided for the outer-loop closure. Then the static gain of the outer loop was increased 
until the altitude mode of motion was neutrally stable. In this way the minimum period of 
the altitude mode of motion that could be achieved with the first-level pilot model was 
determined. 

For aircraft with low values of determining the best response was more 
involved for two reasons. First, the frequency of the angle-of-attack mode of motion 
was low, and the inner-loop closure had to be used to increase this frequency to an 
acceptable value. Second, the outer-loop closure destabilized the angle-of-attack mode 
of motion. Therefore, a compromise was necessary to achieve a suitable system 
response. 

The results of these trial-and-error procedures are shown in figure 4, where the 
periods of the two modes of motion are noted in the plot of as a function of 2£co n 

^in the location of and 2?u> n for which they were obtained^ together with the pilot- 

rating boundary from reference 10. The periods of the modes of motion are given 

instead of the frequencies because it is believed that the periods could be more easily 

o 

visualized and would therefore be more meaningful. For high values of c 0 ^ the period 
of the altitude mode of motion varied from 10 seconds when 2£u> n was low to 4 seconds 
when 2£> n was high. The period of the angle-of-attack mode of motion was very short 
over the entire range of 2£co n , the indication being that control on angle of attack was 
very tight. From these results an altitude -mode -of -motion period of about 5 seconds 
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Figure k.- Results of trial-and-error procedure for obtaining best closed-loop 
system response for altitude control, with first-level pilot model used. 
(Values of and P a are in seconds; the pilot-rating boundary from 

ref. 10 is presented for comparison. ) 


was arbitrarily selected as being a suitable value to use to separate satisfactory from 

9 

acceptable system response. For low values of u>^ it was found that a 5-second 
altitude-mode-of-motion period could be obtained in combination with an angle -of- attack - 
mode-of-motion period of 2 or 3 seconds. It was believed that this angle-of-attack-mode- 
of-motion period was in the range of closed-loop system response that would be desired 
by a human pilot for an inner-loop variable. 

2 

It was also concluded from these calculations that in the plot of cu n as a function 
of 2£co n , there were two distinct regions with different combinations of critical closed- 
loop factors. For the high values of the period of the altitude mode of motion and 
the damping of both the altitude mode of motion and the angle -of -attack mode of motion 
appeared to be the critical factors; the period of the angle -of -attack mode of motion and 
both the period and damping of the control mode of motion seemed to be satisfactory 
without further consideration. For the low values of the periods of the altitude and 

the angle -of -attack modes of motion and the damping of the angle -of -attack mode of 
motion were the critical factors. 


To determine attitude response, the outer loop of the system shown in figure 3 was 
discarded and just the inner loop was used. The procedure followed was similar to that 
described for altitude response. In this case the dominant system modes were a first- 
order pitch-angle response and an oscillatory angle -of -attack mode. For high values 
of the critical factors were the time constant of the pitch-angle mode of motion and 
the damping of the angle -of -attack mode of motion. For low values of the critical 

factors were the period and damping of the angle -of -attack mode of motion. The results 
of the trial-and-error procedure for obtaining the best system response are shown in 
figure 5. For high values of the pitch-angle time constant varied from 6.3 seconds 
when 2£co n was low to 2.6 seconds when 2£a> n was high. When a>2 was low, the 
period of the angle -of -attack mode of motion was around 2.5 seconds. 
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Figure 5»- Results of trial-and- error procedure for obtaining best closed- loop 
system response for attitude control, with first-level pilot model used. 

(Values of T and P a are in seconds; the pilot-rating boundaries from 

refs. 10 and 13 are presented for comparison.) 

Also shown in figure 5 is the pilot-rating boundary from reference 10. However, 
since this boundary was derived from tests in which the task was to control altitude 
(that is, to follow a glide slope), it may not be appropriate for comparing with attitude 
control. Therefore, a boundary published in reference 13 is also presented. The con- 
trol task in reference 13 was primarily long-range tracking at low speed and therefore 
corresponds to the computed conditions. A comparison of these two experimental 
boundaries with the computed results and with results from previous investigations of 
pilot response (which indicated that a human pilot desires a time characteristic of 
around 2 seconds for a single-loop system response) leads to the selection of a pitch- 
angle time constant of 2.6 seconds along with neutral stability for the angle -of -attack 
mode of motion for the boundary conditions for high values of and to the selection 

of an angle -of -attack period of 2.5 seconds and neutral stability for low values of 

With the values of suitable closed-loop system characteristics determined, the 
synthesis method described in reference 14 was used to compute the curves which would 
satisfy these specifications. With this method it is possible to specify any number of 
desired closed-loop system characteristics, provided that an equal number of open-loop 
factors are left free to adjust. Values for the open-loop factors that are required to 
meet the specified closed-loop characteristics are determined. In the present analysis 
the pilot-model static gains are used as free control factors. Also, to define boundaries, 

p 

a series of values for wj: could be given, with the value of 2£ci> n left free to adjust to 

, 11 p p 

a required value. Or 2?co n could be given, and determined. Both w* and 2?co n 
can be related to aircraft stability derivatives, as shown by the following equations: 


,2 _ 


w n "^a^q ~ 


( 12 ) 


2?co n = L a - M q 


(13) 
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Therefore, for altitude control it was possible to specify either of the two following 
combinations of closed-loop system characteristics: 

(1) A period of 5 seconds for the altitude mode of motion and zero damping for both 

the altitude mode of motion and the angle -of -attack mode of motion 

(P h = 5 seconds; ? h =0; = o) 

(2) A period of 5 seconds for the altitude mode of motion and a period of 2.5 seconds 

and zero damping for the angle-of-attack mode of motion = 5 seconds; 

Pq, = 2.5 seconds; = o) 

For attitude control the two sets of closed-loop system characteristics which could 
be specified were the following: 

(1) A pitch-angle time constant of 2.6 seconds and zero damping for the angle-of- 

attack mode of motion (T = 2.6 seconds; = 0) 

(2) A period of 2.5 seconds and zero damping for the angle-of-attack mode of 

motion ^P fl = 2.5 seconds; = 0^ 

It was believed that these combinations of closed-loop characteristics would define 
all pilot-rating boundaries. 

Table in presents some sample solutions for altitude control. Presented is a com- 
plete listing of the open-loop aerodynamic parameters and pilot- model static gains that 
were determined by the computations and of the closed-loop characteristics, both the 
ones specified and the remaining ones that were left free. Solutions are considered 
suitable only if the unspecified closed-loop characteristics are such that they would make 
negligible contributions to the time response of the system. 

TABLE HI 


COMPUTED SYSTEM PARAMETERS FOR ALTITUDE CONTROL 
[Lq, = 0.585 per sec; first-level pilot modelj 


Open-loop parameters 

Closed-loop characteristics 

u 2 

n 

2 C"n 

M a 

M d 

K 0 

K h V 

“h 

2 ?"h 

( o 

a 


“c 

2i >c 

0.63 

6.0 

2.54 

-5.41 

23.8 

2.02 

1.26 

0.43 

2.5 

0 

8.4 

15.4 

1.62 

5.0 

.97 

-4.41 

20.1 

2.14 




36 





8.0 

14.6 

2.62 

4.0 

-.62 

-3.41 

16.4 

2.32 




29 





7.4 

12.7 

3.61 

3.0 

-2.20 

-2.41 

12.6 

2.60 
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6.9 

12.8 

4.62 

2.0 

-3.80 

-1.41 

8.9 

3.15 


f 

.04 

\ 

/ 


f 

6.4 

12.0 

30 

3.2 

-28.5 

-2.62 

38.8 

4.41 

1.26 

0 


5.22 

o 

7.6 

13.2 

25 

3.18 

-23.5 

-2.60 

33.4 

4.31 





4.84 



7.6 

13.2 

20 

3.10 

-18.5 

-2.52 

27.8 

4.19 





4.43 



7.4 

13.1 

15 

2.92 

-13.6 

-2.34 

21.9 

4.02 





3.96 



7.2 

12.9 

10 

2.56 

-8.84 

-1.98 

15.6 

3.78 





3.38 



6.9 

12.4 

5 

1.79 

-4.29 

-1.21 

8.2 

3.38 

\ 

f 

\ 

f 

2.54 

\ 

t 

6.3 

11.7 
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Results 


Effects ot airplane shor t-period characteristics .- The computed curves (determined 
by using the synthesis method) which satisfy the closed-loop system characteristics for 
altitude control are shown in figure 6. The first-level pilot model, as described in the 
section on pilot models, was used. The pilot-rating boundary from reference 10 is also 
presented. 



Figure 6.- Comparison of the curves computed by using first-level pilot model 
(no lead) to provide specified closed-loop system characteristics for alti- 
tude control with the pilot-rating boundary from reference 10. (Values 
of Pft and P a are in seconds.) 

It can be seen that the system requirement Pj^ seconds, £^=0, 1(^=0 defines a 
curve that agrees well with the upper branch of the pilot-rating boundary and that the 
requirement Pj^ seconds, P Q ,=2.5 seconds, ?q,= 0 defines a curve that agrees well 
with the lower branch. The nature of these two curves seems to agree with the detailed 
pilot comments that result from flight tests. When pilots are asked to evaluate an air- 
craft which has a low short-period natural frequency, they complain that the aircraft 
seems sluggish and that they have to overdrive it. This comment is believed to corre- 
spond to the fact that the short period of the closed-loop system is a critical factor with 
this type of aircraft and that the inner-loop closure has to be used to decrease this period 
to 2.5 seconds. When pilots are asked to evaluate an aircraft with a high short-period 
natural frequency, they complain that the aircraft tends to bobble. This comment is 
believed to correspond to the conflicting factors that the inner loop must be closed with a 
high gain to provide adequate compensation for the outer-loop closure but that this high 
gain also drives the inner loop toward instability. These agreements between the shape 
and nature of the computed and experimental curves lead to the conclusion that they both 
describe the same phenomenon. 

Next, a 1-second lead time constant was added to the inner-loop pilot model to pro- 
duce the second-level pilot model, and the synthesis method was used again to define a 

O 

curve in a plot of as a function of 2£c u n . With the addition of lead the critical 

closed-loop factors were = 5 seconds, Pq, = 2.5 seconds, and % a = 0. The curve 
which satisfies these conditions is shown in figure 7, and this curve is seen to agree well 
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with the pilot-rating boundary between acceptable and unacceptable aircraft. From this 
result it is concluded that the use of the second-level pilot model to provide the specified 
closed-loop characteristics can separate acceptable from unacceptable aircraft. 



Figure 7-- Comparison of the curves computed hy using first-level pilot model 
(no lead) and second-level pilot model (lead) to provide specified closed- 
loop system characteristics for altitude control with the pilot-rating 
■boundaries from reference 10. (Values of P h and P a are in seconds. ) 

The boundary defined by the specification Ph=5 seconds, ?h=0, £ a =0 was also 

computed with pilot lead and was found to lie in the same area as the boundary described 
in the previous paragraph. However, the solutions did not meet all the requirements for 
this specification because the period of the angle -of -attack mode of motion, which was 
not specified, was longer than 2.5 seconds. Therefore this boundary was not plotted. 

With the linear, constant-coefficient pilot model used in these computations, there 
is no restriction in the region of low and negative and therefore the cutoff for 
acceptable aircraft near zero frequency that occurs in the experimental data is not pre- 
dicted. It would appear that some consideration of leaving the aircraft uncontrolled for 
some length of time should be included in predicting pilot ratings. A consideration of 

O 

this type would provide a boundary near u> “ = 0. 

Next, attitude control determined by using just the inner loop, with no pilot lead, 
was examined. As can be seen in figure 8, the closed-loop specifications of 
T=2.6 seconds, and P Q ,=2.5 seconds, ? Q! =0 result in two curves which, again, 



Figure 8.- Comparison of the curves computed by using first-level pilot model 
(no lead) to provide specified closed-loop system characteristics for atti- 
tude control with the pilot-rating boundaries from references 10 and 1J. 
(Values of T and P a are in seconds. ) 
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agree well with the pilot-rating boundaries from references 10 and 13. It appears that 
these specifications, together with the first-level pilot model, will also separate satis- 
factory from acceptable aircraft. 

When lead was added to the pilot model for the attitude control task, the surprising 

o 

result was that system response was degraded in the region of high co*, as is shown in 
figure 9. As was true for altitude control, the addition of lead eliminated all restriction 



Figure 9.- Effect of lead in pilot model for the attitude control task. 

(Values of T are in seconds.) 

in the region of low and negative cu*. However, instead of moving the boundary for satis- 
factory control in the region of high to the left, as was the result with altitude con- 

trol, the boundary was moved to the right. This would seem to be a very confusing situa- 
tion in which the addition of lead would help the pilot bring the aircraft under better 
control for altitude control but would degrade the system for attitude control. Pilot 
comments indicate that an exceptional amount of confusion does exist for aircraft in this 
region. In reference 15 pilot comments for an aircraft with a normalized lift-force 
derivative of 0.5 per sec, an undamped natural frequency of 0.7 cps, and a damping ratio 
of 0.3 (see data point in fig. 10) indicate that altitude control was fair whereas attitude 
control was poor because of a tendency for the aircraft to bobble. The computed results 
are believed to be in good agreement with the detailed pilot comments, and this agree- 
ment emphasizes the fact that, for aircraft in the acceptable region, the task has a 
noticeable effect on pilot rating. 



Figure 10.- Comparison of altitude and attitude control 'boundaries 
computed by using second- level pilot models. 
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Effect of T fl o or L .- Another set of results which can be used to show the 

U , Ci Cz 

validity of the pilot model is the data that relate changes in pilot ratings to changes 
in L ft . The effect of on pilot ratings has been summarized by the statement that 

pilots prefer that 

T 0 , 2 w n “ r “ 4 ra <l 

o 

This equation indicates that for = 0.585 per sec, the preferred value for w* is 
o o n 

( T 2 id\ 9 /rad\^ 

) . In figures 6 and 8 it can be seen that co„ = 5.5( | corresponds to that 

sec/ n \secj 

part of the computed boundary that gives the lowest value of aircraft damping 2£cu n for 

satisfactory system response. With ~L a = 1.3 per sec, the handling -quality data indicate 

p /rad 

that the preferred value for c o* increases to 27 ( . The boundaries for the same 

n \,sec/ 

system characteristics, with no lead in the pilot model, were also calculated with 
increased to 1.3 per sec. These results for altitude and attitude control are shown, 
together with the results for = 0.585 per sec, in figures 11 and 12, respectively. 



Figure 11.- Comparison of altitude-control boundaries computed by using 
first-level pilot model for two values of 1^. (Values of P^ and 

P a are in seconds; 10,2°^ is in radians.) 



Figure 12.- Comparison of attitude- control boundaries computed by using 
first-level pilot model for two values of 1^. (Values of T and 
P a are in seconds. ) 
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The compuated system characteristics indicate that the satisfactory area in the plot of 
as a function of 2£cn n has been greatly expanded in the region of high with 
the location of 0 w n = 4 being within this expanded area. The satisfactory region 
is slightly restricted in the region of low o>“ for the higher value of L a . Thus it can 
be seen that the computations are compatible with the trends shown in the experimental 
pilot-rating data. 

With pilot lead added, the results follow the same trends as noted previously. The 

computed boundary for altitude control is shown in figure 13, where the addition of lead 

is seen to move the curve further to the left. A part of this curve is shown dotted for the 

o 

following reason. For high values of cu n the solutions obtained with the synthesis 
method were suitable in that the period of the angle-of-attack mode of motion, which was 
not specified for this boundary, was shorter than 2.5 seconds. However, for low values 

o 

of the period of the angle-of-attack mode of motion became longer than 2.5 seconds, 

and therefore the solutions were not considered suitable and are not shown. However, 
approximate analog solutions obtained by a trial-and-error method and judged on the 
basis of the computed time histories indicated that suitable system response could be 
obtained for aircraft characteristics to the right of the dotted line. These solutions were 
therefore used to complete the boundary. 



2t> n , rad/sec 


Figure 13. - Comparison of altitude-control boundaries computed by using 
first-level pilot model (no lead) and second- level pilot model (lead). 

(Values of P^ and P a are in seconds.) 

With pilot lead the boundary computed by the synthesis method for the specification 
Pk=5 seconds, P a =2.5 seconds, ? Q! =0 also fell in the same area as the boundary which 
specified Ph=5 seconds, £^=0, However, the control mode of motion that 

resulted had a very long time characteristic, and therefore the solutions were considered 
not suitable. 

The computed boundary for attitude control, with L a =1.3 per sec and with 
pilot lead, is shown in figure 14. The results follow the same trend as noted for 
L^, = 0.585 per sec when pilot lead was added. The addition of lead removes the 
restrictions for low values of but moves the boundary to the right for high values 
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of ^n- However, the boundaries are less restrictive for L a = 1.3 per sec (fig. 14) 
than they were for Lq, = 0.585 per sec (fig. 9). 



Figure lit-.- Comparison of attitude- control 'boundaries computed by using 
first-level pilot model (no lead) and second- level pilot model (lead). 

(Values of H and Pa are in seconds. ) 

Flight tests have been used to define pilot ratings as functions of and 2£o; n 
and as functions of T 0 2 and a> n . The previous paragraphs have shown that computed 
boundaries agree well with these flight data. The combined effects of cu n , 2t> n , 
and Tg 2 have not been systematically studied in flight tests. One of the immediate 
advantages to be gained from the present investigation is that it indicates the effect of 
the combination of all three parameters. In fact, the present method can provide 
answers to the question of pilot ratings for any number of system parameters and can 
thereby alleviate the need for extensive flight-test programs. 

Effect of control -actuator dynamics .- A third influencing factor on pilot ratings that 
can be checked with the use of pilot models is the degrading effect of inserting control- 
actuator dynamics in the control system. It has been shown that if the actuator has a 
natural frequency of less than 20 rad/sec, the pilot ratings are sometimes drastically 
affected. Therefore, system characteristics were computed with a second-order actuator 
system inserted in the pilot-aircraft system. The actuator was assumed to have a 
natural frequency of 10 rad/sec and a damping ratio of 0.7, which is the type of actu- 
ator response that is very likely to be encountered in a power-boost control system. 
Attitude control with no pilot lead was examined. For the low value of L ff (that is, 

Lq, = 0.585 per sec), the boundaries computed with no actuator (perfect response assumed) 
and with the 10-rad/sec actuator are presented in figure 15. It can be seen that the addi- 
tion of the actuator restricts the satisfactory region, or conversely, for given aircraft 
short-period characteristics the pilot rating would be degraded. In reference 12 two low- 
speed aircraft configurations (configurations A and B) were tested with variations in the 
actuator natural frequency, and the locations of these aircraft characteristics are also 
plotted in figure 15. For configuration B with a very high actuator frequency (essentially 
perfect response), a pilot rating of 2 (satisfactory) was given. For the same aircraft 
configuration combined with an actuator frequency of 10 rad/sec, a rating of 3 or 3.5 
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Figure 15. - Effect of actuator dynamics on attitude- control boundaries computed 
by using first- level pilot model (no lead) with = O.585 per sec. (Values 
of T and P a are in seconds.) 

(borderline rating between satisfactory and acceptable) was indicated. The computations 
show that this aircraft is inside the satisfactory region (a rating of better than 3.5) with 
perfect actuator dynamics and moves outside the satisfactory region when the actuator 
with the frequency of 10 rad/sec is added. The experimental and computed results are 
therefore seen to be in good agreement. For configuration A the pilot ratings changed 
from 5 to 6 with the addition of the actuator, and this change again agrees with the com- 
puted results, although the overall level of the experimentally derived ratings is more 
towards unsatisfactory than the computed results indicate. It is pointed out in refer- 
ence 12 that some extenuating circumstances were present that were believed to degrade 
the overall pilot ratings. 

For the high-speed condition (L a = 1.3 per sec), the computed boundaries are pre- 
sented in figure 16. Again, changes occur in the position of the boundaries when the 
actuator is added. The locations of the characteristics of two high-speed configurations 



Figure l 6 .- Effect of actuator dynamics on attitude- control boundaries computed 
by using first-level pilot model (no lead) with 1 ^ = 1.3 per sec. 

tested in reference 12 (configurations E and F) are also presented. For configuration E 
the ratings presented in reference 12 changed from 2 to 5 with the addition of the actuator. 
The computations show that the configuration is well inside the satisfactory region with 
no actuator and just on the boundary with the 10-rad/ sec actuator, so fairly good agree- 
ment exists between the computed and experimental results. For configuration F no 


23 


change was noted in the experimental test results, and the computations show that this 
configuration falls in a region where the addition of the actuator has very little effect. 
Again, there is good agreement between the experimental and computed results. 

APPLICATIONS 

Since the pilot-model levels and the desired closed-loop system characteristics 
defined in the preceding sections of the paper have been shown to provide a good corre- 
lation with a very wide variety of experimentally determined pilot ratings, it is con- 
cluded that these pilot-model and closed-loop system criteria can be used to predict the 
pilot rating of the longitudinal handling qualities of proposed aircraft. The rating is pre- 
dicted either by determining the values of the lead and lag time constants required in the 
pilot model to achieve the specified closed-loop characteristics or by determining if these 
closed-loop characteristics are achieved with a given pilot model. 

It is believed that the same general procedure could be followed to determine the 
lateral handling qualities of an aircraft or the handling qualities of other types of vehicles. 
However, it might be necessary to adjust the specified closed-loop characteristics to 
correspond to the level of system response to which human pilots have become accus- 
tomed. This adjustment, of course, could be made by determining the best correlation 
with past handling-quality data for the particular type of vehicle under consideration. 

CONCLUDING REMARKS 

System response characteristics obtained by using pilot models have been shown to 
correlate with pilot ratings as follows: 

(1) In the shape of the boundaries in the plots of airplane short-period 

characteristics 

(2) In the change in pilot ratings with change in normalized lift-force derivative 

(3) In the change in pilot ratings with the insertion of control-actuator dynamics 

It is therefore concluded that a prediction of pilot ratings can be made by determining the 
level of response of the pilot required to achieve specified system characteristics. This 
method can also be used to compare proposed airplane designs and to provide a rational 
engineering interpretation of pilot comments. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., December 7, 1970. 
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APPENDIX A 


SYSTEM PERFORMANCE DETERMINED WITH 
PILOT TRANSFER FUNCTIONS 

Because of the great interest in predicting system performance, an example of the 
response of a pilot-model — aircraft system to a glide-slope command is presented. 

This computed time history is compared with a time history obtained from an actual 
flight test. 

The flight-test data were taken from reference 10 and are shown in figure 17. 

Figure 18 was computed by using the second-level multiloop pilot model, as described in 
this paper, and aircraft characteristics approximately the same as those in the flight- 
test data of figure 17. The pilot model used also included a remnant added to the pilot- 
model output. This remnant was obtained by taking a white-noise source, passing this 
white noise through a filter that was identical to the response characteristic of the pilot, 
and adding it to the output of the pilot model. The expression for the remnant is there- 
fore K n /(1 + 0.2s)^. The amplitude of the remnant K n was adjusted so that the mean- 
square value of the remnant was approximately 0.5 of the mean-square value of the total 
pilot-model output. It has been shown in reference 2 that a pilot remnant defined in this 
manner is appropriate. 

The time history in figure 18 can be seen to agree well with the time history 
obtained from the flight tests (fig. 17). The general amplitude and frequency of the angle- 
of-attack and pitch-angle time histories are in good agreement. One difference occurs 
in that the computed results show a large response at the discontinuous point in the com- 
manded flight path, which the human pilot was able to anticipate and smooth out. 



0 10 20 30 40 50 60 

Tlme,f, sec 

Figure 17 . - Records obtained from reference 10 for flight test of 
a landing approach, = 2.30 rad/sec; £ = 0.20. 
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APPENDIX A - Concluded 



Figure l8.- Landing- approach time history com- 
puted by using second-level multiloop pilot 
model with remnant, = 2.5 rad/sec; 

£ = 0.32. (To facilitate comparison with 
fig. 17^ Ax and A0 are given in degrees 
as well as radians.) 




Time.t, sec 


Figure 19- - Landing- approach time his- 
tory computed by using second- level 
multiloop pilot model without 
remnant, cdq =2.5 rad/sec; £ = 0.32. 
(To facilitate comparison with 
fig. 17, Ax and A0 are given in 
degrees as well as radians.) 


To illustrate the effect of the remnant on the system response, figure 19 was com- 
puted without the remnant. The longer period, lightly damped altitude mode of motion 
is more apparent in figure 19 than it is in figure 18. 
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